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Abstract 

Some techniques for the study of intermittency by means of wavelet transforms, are presented 
on an example of synthetic turbulent signal. Several features of the turbulent field, that cannot 
be probed looking at standard structure function scaling, become accessible in this way. The 
concept of a directly measurable intermittency scale, distinct from the scale of the fluctuations, 
is introduced. A method for optimizing the analyzing wavelets, which exploits this concept, and 
allows to minimize non-local contributions in scale to wavelet correlations, is described. The 
transition from a wavelet to a Fourier transform based description of an intermittent random field, 
and the possibility of using Fourier correlations to measure intermittency are discussed. Important 
limitations in the ability of structure functions to give a local in scale description of intermittent 
random fields, are observed. 
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I. Introduction 

Structure functions and scaling arguments, are basic tools in the study of turbulence. One of 
the key issues in the field, intermittency, is expressed in this language through the statement, that 
the scaling of the structure functions Si{q) ~ {\Aiv\'^), with Aiv{x) — v{x + l) —v{x), is non-trivial. 
This means simply, that in high Reynolds number turbulence, there is a large "inertial" range of 
scales where Si{q) obeys a power law: Sq{l) ~ Z^', with (q a nonlinear function of the order q. 

Many explanations have been proposed for this phenomenon, but to date, nothing definitive is 
still available (see for a review). There are several reasons for this state of affairs; one is perhaps, 
that scaling by itself is not able to provide a sufficiently complete description, of what is happening 
in the turbulent field. Within Kolmogorov 1941 theory |^, energy conservation is able by itself to 
fix a value of the scaling exponent: Q"' — 9/3, which is an acceptable lowest order approximation 
for Cg. However, no comparable symmetry based argument has been proposed, which was able to 
fix the value of the intermittency correction C^"' — Cq. 

The fact is that, although the mechanism of intermittency generation is likely to be universal, 
possibly associated with some property of the nonlinear energy transfer in the inertial range, it 
is also likely that its effect is not exhausted in the production of scaling corrections. In fact, 
very different kinds of intermittency, like e.g. the one associated with the long and thin vortices 
observed in numerical simulations , and the one that would be obtained in a random beta-model 
picture of turbulence B, produce equally acceptable spectra of scaling corrections. 

Wavelet analysis |^ has often been proposed as an alternative tool in the study of turbulent 
intermittency. The possibility of having an additional degree of freedom associated with the wavelet 
shape, beyond position in space and scale, makes these objects, particularly appropriate to study 
coherent structures and the geometrical properties of intermittency However, the great freedom 
to describe geometrical features in two- or three-dimensional settings, thus becoming available, has 
hindered perhaps an extensive use of these techniques in turbulence theory. 

An alternative use of wavelets has been to employ them as building blocks in the generation of 
artificial turbulent signals, trying to reproduce the kind of velocity time series one gets in experi- 
ments. This kind of technique, first introduced in Q to study connections between dissipation and 
velocity intermittency, and inH , has been later developed systematically in |^ (see Q for recent 
references). More recently |11|], this approach has been used to provide a kinematic explanation 
for the kind of energy spectra developing in wall turbulence. The picture of a superposition of 
eddies at different scales obtained in this way, appears particularly natural to study intermittency 
effects. However, even in such one-dimensional settings, there is a great freedom in the choice of 
the eddy generation mechanism. It is clear, for example, that the same spectrum of anomalous 
exponents can be obtained, both from a "random eddy model" with inclusion of intermittency, 
like the one considered in ||ll| , and from a multiplicative cascade of the kind described in |^ . (To 
get such an identically scaling signal, out of the first model, it is enough to randomly permute the 
wavelets at the different scales in the second). Now, different mechanisms of turbulence synthesis 
imply, to some degree, different assumptions on the real turbulent dynamics. It is therefore of 
some relevance, to devise methods which allow to identify these mechanisms, from the statistical 
properties of the signal. 

In light of the present discussion, synthetic turbulence appears to be the appropriate "training 
facility" , in which to test different techniques of wavelet analysis, and their ability to reveal specific 
intermittency feature s ||T^ . 

Recently, spatial |13[ and scale correlations between wavelet components, have been used 
to probe the cascade structure of turbulent signals. Here, the interest is focused on two different 
issues: the choice of the wavelet in both the signal generation and analysis, and the amount of 
phase space available to each "building block" wavelet in the generation algorithm. This phase 
space is given by the relative position and scale of the generated wavelet, with respect to the parent 
one, and in typical algorithms of signal generation it consists just of a single point. It turns 
out that both issues of wavelet shape and phase space availability, have important consequences 
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as regards the ability of a structure function to detect features of the turbulent dynamics that are 
really local in scale. 

In the next section, the generation algorithm for the synthetic turbulent signal is introduced, 
and in section III, the basic wavelet structure function properties are derived. A systematic analysis 
of the turbulent statistics dependence, on the shape of both the analyzing and the building block 
wavelets, is carried on in section IV. Given the one-dimensional nature of the signal, geometrical 
aspects are minimal, and this leaves out only one essential degree of freedom in the choice of the 
wavelets: their "number of wiggles", i.e. the product of their dominant wavevector and spatial 
extension. Section V contains discussion of the results and conclusions. We leave in the appendix, 
the analysis of the case in which the cascade is probabilistic in space and discrete in scale. 



II. Synthesis of an artificial turbulent signal 

Introduce the Gaussian wavepacket: 

ws(k,y,x) ^ exp{ik{x - y) - {x - y)'^/Xs); Xs^as/k (1) 

and let: 

= XlXl^h,.u;s(fch„,yh„,a;) (2) 

n h„ 

be a real random field, which should mimic the time signal from a fixed position velocity mea- 
surement in a turbulent flow. Following standard practice the building block wavelets ws are 
generated through a cascade process, to model the mechanism of energy transfer in the turbulent 
flow. The vector index h„ = {Hq, hi, ....hn) identifles then the position of the wavelet in the cascade 
through the sequence of its ancestors: the integer /i„ labels the /i„-th daughter wavelet generated 
at the n-th step in the cascade, by wavelet h„_i. 

Notice, however, that an intermittent random field could be generated, without any reference 
to cascade processes, either by varying appropriately the space density of the wavelets with scale, 
or by making the distribution of the amplitudes A more intermittent as k grows, but keeping the 
wavelets randomly distributed in space . 

From reality of for each index h„ with positive components, there is a wavelet with index 
— h„ such that A_h„ = A^^^ yh„ = y-h„ and fch„ = — fc-h„- The cascade is assumed to be local in 
n in the sense that the probability that a given wavelet h„ has a certain value of its parameters 
^h„ = {In Ah„ , In fch„ , 2/h„ }, can be written in terms of transition probabilities, as: 

P{^h, ^ ?hj = j da,-deh„_,p(ao ^ ax)-P(eh„_i ^ aj. (3) 

For the sake of simplicity, the transition probabilities are assumed to factorize into their In A, In k 
and y components, with scale invariance forcing the cascade to be governed by a multiplicative 
random process: 

P(C ^ O = PA{A'/A\k'/k)pk{k'/k)p,{k\y ~ y'\). (4) 
The relative phase of A' and A is assumed random, and we take: 

J d\n ypA{y\x)yP = CpX''^", (5) 

in order to get power law scaling in the structure functions. Given Eqns. (4-5), the transition 
probability over n steps is in the form: 

P{£,h^ £,hp) ^PA{p~n,AhJAhJkhJkhjPk{p-n,khJkhjP^{khJyh^ -2/h„|)- (6) 

At each step n in the cascade, the wavelets distribution in scale is peaked at fc„ = fco exp(nz), with 
the wavevectors kg distributed around a characteristic large scale L: (ko) = L^^. Each mother 
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wavelet generates exp(^) daughters; this insures that the mean degree of overlap between wavelets 
in k-y space is scale invariant. 

There are several reasons to consider a mechanism of turbulence synthesis, in which the wavelets 
are distributed in space and scale in a probabilistic way, rather than on a rigid lattice. The main, 
rather "philosophical" motivation, however, is to try considering the building block wavelets, more 
like eddies (or components of bunches of eddies if as is large) , than like basis functions at fixed 
position in space; this also in view of possible extensions of the model to time dependent situations, 
in which the eddies are mobile. 



III. Analysis of an artificial turbulent signal 

The choice of analyzing wavelet is in general arbitrary. A Gaussian wavepackets, however, has 
the minimum spread in fc — y space and allows to retain the maximum simultaneous information 
possible in space and scale. We thus take from the start the analyzing wavelets to be derivatives of 
Gaussian wavepackets: {—k)~^d^WA{k,y,x), for which in general: = A/c 7^ as- The components 
of ^'(a;) with respect to this set of wavepackets are defined as follows: 

■^kv = X-Ak-' j dxw*A{k,y,x)dl^{x). (7) 

together with the associated structure functions d^'fej/l')- In order for these structure func- 
tions not to be dominated by the largest scales in '^{x), given standard Kolmogorov scaling 
for ^, it is necessary that the parameter s in Eqn. (7) be at least equal to one. To calcu- 
late them , we need to evaluate first the component of a building block wavelet on an analyzing 
wavelet. Because of random phase of A, we will need only its square modulus: C{k,y;k' ,y') = 
l^-m J w*^{k,y, x)d'^Wg{k' ,y' , x)\^ . If the probabilities Pk vary sufficiently slow, and s is large 
enough to kill the contribution to the A;-integrals from k small, it is possible to write from the start: 

C(.„ . exp ( - (A.' + -||#^!) ) (8, 

with Ak = kA — ks and Ay = yA — ysi a-nd the integrals over k, that arise in the averages involved 
in the correlations can be carried out by steepest descent, disregarding the contribution at ~ L~^. 

The simplest correlations are {\^ky\^) and (|^'fej;p|\E'fe'y'p). Given the random phase of A, it is 
easy to see from Eqns. (2) and (8), that a 2n-ordcr correlation will receive contribution at most 
by n eddies, and the 4-th order correlation will be in the form: 

= ^ F„(i^h„ 1^) + E G«p(i^h„ n ?). (9) 

n np 

Prom the form of we have power laws for the amphtude correlations; 

(|AhJ')=C2(A;h„L)-^^ (10) 

and 

{\A^S\Ai.'S) = C4(fch;.L)-'^*(fch„/fch;.)''^Mfch;.AhJ^^-'^= (11) 

with A;h„ > kii'^ and p the cascade step at which the genealogical tree of h„ and branches: 

hi = h'^ioT Q <i < p and hi 7^ h[ for i > p. Thus, the lower the branching takes place in the tree, 
the closer the correlation gets to its disconnected limit (|Ah„ P|Ah;„ |^) = C4{kh'^kh^L'^)~''^ . 
Prom Eqns. (2) and (7) we obtain, for the second order correlation: 

(|*fcj/|'> =C2 / dyd\nkY,kn{Pk{n,k/ko)){kL)-^'C{k,y;k,y), (12) 
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where {Pk{n,k/ko)) is averaged over ko (from summing over ho in h„), and C{k,y;k' ,y') the 
square modulus of the component of wavelet ws{k\y' ,x) with respect to the analyzing wavelet 
WA{k,y,x): C{k,y;k',y') = \k-"'Jdxw*A{k,y,x)d^w*sik',y',x)\^. 

We have an analogous expression for the one-eddy contribution to (|^'fcj,|^|^fe/y'p): 

(l*fe!/l'l*fcVl')l 

= C4 f dydlnfc^fc„(Pfe(n,fc/fco))(fci)-^*C(A;,2/;fc,y)C'(/e',y';fc,y). (13) 

n 

For the two-eddy contribution, we have instead, for k' > k: 

/An k 
dlnMlnfc'dydy' / dlnk{kL)-^''{k'/k)-^^{k/k)^*-'^^^ 

^ E E E ^{PkiP, k/ko))Pk{n - p, k/k)Pk{m - p, k'/k) 

n m=0 p=0 P 

xP^{k\y-y'\)C{k,y-k,y)C{k',y'-~k',y'). (14) 

where iJim/kp)Px(k\y — y'\) gives the space density at y' of eddies generated from the branching 
at hp, given the presence of an eddy h„ at y. The factor km/kp = exp((m — p)z) is the actual 
number of eddies generated from the branching at hp. The factors kn entering Eqns. (12-14), 
conversely, give the space density of wavelets of typical size as/kn, at the n-th step in the cascade. 

The cascade structure is characterized by a discrete component through the sums entering 
Eqns. (12-14). Since the cascade steps are independent, the width of the cumulative distribution 
Pk{n,k'/k) is niAz with Az the width of pfc: Az'^ = J dlnfc'| Ink'/k - z\^pk{k' /k). Thus, if the 
separation of the scales entering the Pk involved in the sums in Eqns. (12-14) is large enough, the 
effect of discreteness will be negligible. The same will occur if Az/z itself, is large enough. In the 
other limit, when Az/z and n are small, oscillation with period z in In k'/k and InkL (lacunarity) 
are to be expected in the scale dependence of the correlations. 



IV. Structure functions and optimization of the analyzing wavelets 

The length Ag introduced in Eqn. (1) plays an important role in the generation of an intermit- 
tent random field, since it identifies the intermittency scale of fiuctuations of size k~^. If as is large, 
indeed, it is not the single fluctuation that is intermittent, but the amplitude of a whole bunch of 
them, extended over a length As. As a consequence of this, one expects that, if As approaches the 
size of the domain, intermittency is lost and a standard, random phase generated gaussian field is 
obtained. Conversely, if one carries on the same operation with the length A^, the final result is 
that, instead of dealing with structure functions, one ends up working with correlations between 
Fourier components of the random field. Again, intermittency is expected to be lost. If one is 
interested in studying intermittency by structure function scaling, the optimal choice should be 
therefore: A^ ~ As- For this reason, it becomes necessary to study the dependence of the structure 
functions {\'^ky\^) on the parameters ua and as- 

If Az/z is not too small, discreteness effects in scale can be neglected and the sums in Eqns. (12- 
14) can be approximated by integrals. The transition probability Pk is in the form: Pk{n,k' /k) = 
/(n. In fc'/fc — nz) ~ f{z~^ Ink'/k, Ink' /k — nz). Hence: 

y^Ckn/k)Pkin,k/ko) ^ / dxz-^e^f{z-Hnk/ko,x) = -(^ , (15) 

where e = lim„^co ln(exp(x)|n) — 0{z~^Az'^) and {■■■\n) indicates average over f(n.x). We 
thus see that scale uncertainty in the process of eddy generation contributes to scaling in such a 
way that Cg ^ Cg = Cg + £• 
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At this point, all the terms in Pk and associated sums drop off Eqns. (12-14) and it is possible 
to evaluate the integrals for k = k' and y = y'. From Eqns. (12-13), we get immediately: 



C2TTas , 



{kL)-i^ and {\^ky\^)i 



C47rag 



(fci)-^i 



A) 



(16) 



Performing some power counting on Eqn. (14), however, wc discover that the dfc integral is 
dominated more and more by large scales the closer we are to trivial scaling: ^4 = 2(2- In this 
regime, ii Px is sufficiently well behaved, it will be possible to approximate Px, from central limit 
theorem arguments, with a Gaussian: 



Px{k\x-x'\) 



k-\ 



,.'|2 



TT^bas 



exp I 



(17) 



with the parameter b characterizing the spatial non-locality of the cascade. Substituting into Eqn. 
(14), we obtain: 

{\^ky\') = {\^kyW + {\^ky\')2 



C47ra| 



(kL) 



< 



asaA 



a. 



27r2 
z^baf. 



In kL 



dx{ 1 



-2x 



b^al 



.(Cl-2C2)^ 



(18) 



This equation is our main result, and tells us how the various eddies in the synthetic turbulent 
field, contribute to structure function scaling. 

The second term in square brakets in Eqn. (18) is the two-eddy part of (l^'feyl^), giving 
for each x = \nk/k, the logarithmic distance of the common ancestor of the two eddies from 
the scale k. For ^4 — 2(^2 small, which is true in the case of turbulence, this integral receives 



contributions from max(0, ^ In "^2^°^ ) < x < InkL, where the integrand can be approximated by 

exp((C4 — 2(^2)^)- We have then the important result, that for kL < exp((2^2 ~ CD^^); which is 
a very large range of scales, the structure function contains a logarithmic two-eddy contribution, 
which comes right from the largest scales in the random field. We obtain then, for the kurtosis 
K4k,as,aA) = {\^ky\^)-^{\^ky\^y. 



Ki{k,as-,aA) 



2C2-C 



asaA , 27r2i?2C2 C 
—2 2" "I ~ RkL 



(19) 



b d 1 

where: R = max(l, (txt")^)- Only for kL 3> exp((2C2 ~ Ci)~^) we reach pure power law scaling: 

27r3ii2C^-C^ 



K4{k,as,aA) ^ ^{kL) 



2C2-C4 



asaA 



+ 



a% + al z-^has{2C:2 - CD 



(20) 



If, instead of looking at the scaling of Ki^{k, as, a^), we study the dependence of this quantity on 
ttA for k fixed, we notice the presence of a maximum at qa = as- This corresponds to the maximum 
possible overlap between building block and analyzing wavelets; for > ag, an analyizing wavelet 
will feel the effect of many eddies at different position in space, while, for < as, these will be 
distributed at different scales. It is important to notice, as it is clear from Eqn. (19), that this 
effect will be felt also in the measured scaling exponents, that, because of the logarithm, will be 
dependent on ua- This effect will be minimum only at ua = as, when the local in scale, one-eddy 
contribution to the structure function is maximum. Conversely, the importance of the two-eddy 
contribution goes to zero when as is large. 

It is important to stress the importance of the smallness of 2^2 — Ci a^nd of the cascade structure 
of the random field. This causes the slow decay of the two-eddy contribution as gets large. In 



6 



a random eddy model of the kind considered in ||Tl|, the two-eddy contribution would scale like 
(fcL)~^^2 whatever the value of a^, the reason being the lack of correlations among eddies. 

When oa/os becomes large enough, the one-eddy contribution can be disregarded, and the 
kurtosis K4 begins to scale in R, i.e. in the ratio Xs/^a, so that we have, to leading order in k: 

K4k, as, aA) ~ (L/Aa)"^^"'^^ (21) 

In fact, for large, wa probes eddies coming from common ancestors, which can be very distant 
in scale from fc, and which become uncorrelated when > L. From Eqn. (21), (|^fcy|^) obeys 
trivial scaling: {\^ky\'^) ^ (fcL)~^'»2 and the intermittent nature of the random field is lost. 

In this regime, the ratio of the one- to two-eddy contribution to (l^fcj^l'') scales, again to leading 
order in fc: ^ 

The factor (kL)^^ in this equation has a very important interpretation when we consider that, for 
fixed Xa ~ L, and kL large, XA'i'ky ^fc, which is the Fourier transform of '^{x) in a box of size 
L. We have then, that Eqns. (21-22) take the form: {\^k\'^) - aLk-^{Lk)-'^^ + hL'^k-'^{Lk)-'^^'' 
which is nothing else than the expression {'^ ki'^ k^) — S{ki+k2 + k3 + k4)C4{ki,k2,k3,k4) + 

6{ki + k2)5{k3 + ki)C2{ki)C2{k^) for ki = -fcz = k^ = -ki with ,5(0) ~ L. Thus, {\^ky\^)2 kill 
(l^feyl^)i by a factor kL that is the term k5{0) that gives the ratio of the disconnected to the 
connected Fourier 4-point correlation of a uniform random field. All this suggests that, perhaps, 
structure functions are not the most appropriate objects to probe features of the field which are 
local in scale; rather, the multiscale Fourier correlation: (^fc^^'fe^^'/cg^'/c^) with ki ^ —kj V i,j and 
ki + ... + k^ = should be taken into consideration: 

{^S,^...^,^)^27T6{k, + ... + k,)^ I ^(fcL)-';^exp(-^^(fc„-fc)2). (23) 

n—l 

Only in this way, it would be possible to avoid logarithmic corrections in the scaling of correlation 
functions, coming from two-eddy effects. 



V. Conclusions 

The motivation for the interest in synthetic turbulence has often been, more in the " output" , 
i.e. in the turbulent field or turbulent signal being produced, than in the dynamical meaning of 
the adopted algorithm. A typical application has been, for instance, the possibility of controlling 



velocity spectra, in the study of turbulent diffusion |15 . When it comes to an issue like intermit- 
tency, however, the interest is more in the algorithm itself and in the effect that different choices 
for it, would have on the turbulent statistics. In the present research, the main result is the ease, 
with which some innocent looking algorithms for the generation of synthetic turbulence, lead to 
strongly non-local effects in structure function scaling. This, despite the local nature of the cascade 
mechanism on which the algorithms are based. 

The only way this phenomenon can be explained is through the interaction of the cascade nature 
of the algorithm, with the probabilistic distribution in space and scale of the eddies. In particular, 
had the eddies been distributed at random, without a cascade structure, the two eddy contribution 
to a 4-th order structure function, would have been trivially: J dkdk'dydy'C{k, y; k, y)C{k, y; k' , y') 
X (fcfc')~'»*'/^ ~ with C{k,y]k,y) the square wavelet component of an eddy [see Eqn. (9)]. 

On the other hand, if the cascade had been rigid, organized on a lattice structure as in this 
contribution would have been simply NC{k,y;k,y)'^k~'''^ ~ fc^''*, with N the number of wavelets 
ws in the lattice, overlapping with the wavelet wa, and ky the typical scale and coordinate of the 
wavelets ws. 
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All these effects manifest themselves in the behavior of structure functions, through logarithmic 
scaling corrections carrying information on the largest scales of the signal. One result of the present 
study could therefore be that, if one wanted simply to generate a random field with prescribed 
multifractal statistics, it would not be a good idea to include probabilistic effects in the space 
and scale distribution of eddies. If, on the other hand, one thinks that all this may have some 
relevance for real turbulence, the conclusion is that there exists a kinematic effect, which could 
limit the ability of objects like structure functions, to detect scale-local features of the turbulent 
field. This confirms the current thinking on the subject, which leads to expect non-locality, among 
the other things, from the convolution nature, in Fourier space, of correlations of order greater 
than two. It is worth stressing the kinematic nature of our result, which is due to the way in 
which contributions from different eddies sum up in structure functions: the dynamics of the eddy 
generation mechanism remains strictly local in scale. 

The question at this point is how to verify that the non-locality effects we have discussed are of 
any relevance in high Reynolds numbers turbulence. An experimental test could be the dependence 
of the wavelet structure function scaling exponents, on the parameter a a- the number of wiggles 
of the wavelet. In this case, the suggestion from the present research, is that there should be a 
value of a A, for which non-local effects become minimal, and which is identified by the maximum 
at fixed scale of the generalized kurtosis [see Eqns. (19-20)]. This maximum has a nice physical 
interpretation in terms of resonance between the analyizing wavelet and intermittency, with the 
wavelet dominant wavevector identifying the scale of the turbulent fiuctuations, and the wavelet 
extension giving the lengthscale over which these fluctuations, act coherently to generate intermit- 
tency. If the picture of a probabilistic cascade is right, however, only a description based on Fourier 
correlation scaling, could allow elimination of this kind of non-local effects, once all disconnected 
contributions to the correlation, are eliminated by appropriate choice of the wavevectors. 

The results of this study are of rather general validity, showing that there is a rather broad 
class of intermittent random fields, for which techniques based on the analysis of structure function 
scaling, are of limited use. If the interest is more in the modelling of turbulence intermittency, 
however, our non-locality effects may be considered more as an artifact of the algorithm of tur- 
bulence generation. On the other hand, it is difficult to a priori exclude a probabilistic cascade, 
as opposed to a "rigid" one, and several arguments in favor and against both are easy to find. A 
drawback of a probabilistic cascade is that the eddies can overlap in k — y space, while being treated 
as independent objects in the random multiplicative process. One may argue back, however, that 
we are dealing with a one dimensional section of a three-dimensional turbulent field, and that 
these overlaps are therefore irrelevant. Conversely, it is more aesthetically pleasing to distribute 
the wavelets in the random field, freely in fc — y space, but then one looses the property of the 
wavelets, of being base functions for the random field. 

In any case, there are practical reasons for being interested in probabilistic cascades. One is the 
possibility of studying the effect of space and scale non- locality in the eddy generation , which 
have a direct interpretation in terms of properties of the energy transfer in real turbulence. The 
second is, that the lack of constrains over the eddy position allows, in a time dependent situation, 
to put these "eddies" in motion, accounting for sweep in time correlations, in a much easier way 
than using fixed wavelets. 

Aknowledgements: I would like to thank Jean- Francois Pinton, Sergio Cilibcrto and Jens Eggers 
for interesting and stimulating discussion. 

Appendix: The effect of discreteness in the cascade 

For the sake of completeness, we examine the discrete limit of Eqn. (14), analyzing how, 
restriction of the phase space available to the eddies at their birth, modifies the scale non-local 
character of the struture functions. 
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An argument analogous to the one used to arrive at Eqn. (17), leads to the expression for Pfe: 



n(„,t7*) = -^exp(-»^H^). (^1, 

The continuous approximation, used to arrive from Eqns. (12-14) at Eqns. (16) and (18), can 
still be applied to the sums in n of Eqns. (12-13) and in p of Eqn. (14), also for Az/z small. 
This, provided Az^ In kL/z > 1. The effect of discreteness remains therefore in the sums over 
n and m in Eqn. (14), which, for Az/z small, are dominated hy n = p + lnt{z~^ In k/k) and 
m = p + lnt{z~^ In k'/k). Prom Eqn. (14), we thus obtain the following limit expression for 
(l*fe.l')2: 

dfc {kL)-'^'i{k' /k)-^'^{klkf^^^^'^ 



fc2 (ln(fc/fc)ln()fcVfc))2 



xP,{k\y - y'\)C{k, y; k, y)C{k, y; k' , y') exp ( - + -%^)) (^2) 

V Az^ ^\nk/k Infc'/fc^^ 

where Akk' = Ink' /k — int{z~^ Ink' /k) is the decimal part of z~^ Ink' /k. This equation differs 
from the corresponding formula for the continuous limit, Eqn. (18), because of the exponential term 
in A^g and A^^, . This term has a fast dependence on k, which can be treated by decomposing the 
dfc integration as: J ^ ~ ^Sn-- / '^^hk where the integration limits in J dA^^ are approximated 
by ±00 thanks to the smallness of Az. After the Gaussian integral in dAj^^ is carried out, the 
remaining sum can be approximated back to an integral: ^ / X- "^^^^ means that in Eqn. 

(A2) we can substitute: 

V Az^\\nk/k \nk'lk)> ^\n{k/k)\n{k' /k)^ V Az'^ In kk' /k"^^ 
Substituting into Eqn. (A2) and using Eqn. (17) we obtain: 

asttA 



Zd A 



a% + a\ 



oi pin kL , 2,2 -/2i2\,i 

^-ziAzbasJo ^ b^al ) V^2a\alAz^) 



This expression differs from the continuous limit described by Eqn. (18), because of the factor 
[x + a'^O - If the discrete approximation must work over the whole domain of integration 

— / 2 _|_ 2 ■s 

in X, it is necessary however that 2a" Az^ ^ Infci. (Incidently, this tells us that this limit is of 

scarce practical interest in a probabilistic cascade, since it is not particularly interesting to have 
an error kAz in the cascade step, which is much smaller than the wavelet spectral width k/as)- In 
the Az limit we get then the following expression for the kurtosis, for InRkL < (2^2 ~ Ci)"^- 



Kiik,as,aA)^^{kLf^'^-^'^ 



TTCo 



asaA , 27rasaAi?2C2 C 

■ In RkL 



+ a\ z^bia\ + as)i 
Logarithmic corrections to scaling remain therefore, also in the discrete limit. 



(^5) 



9 



References 

[1] U. Frisch, Turbulence (Cambridge University Press, Cambridge, England 1995) 
[2] A.N. Kolmogorov, C.R. Acad. SCI USSR 30, 299 (1941) 
[3] A. Vincent and M. Meneguzzi, J. Fluid Mech. 225, 1 (1991) 

[4] R. Benzi, G. Paladin, G. Parisi and A. Vulpiani, J. Phys. A (Math. Gen.) 18, 2157 (1985) 

[5] I. Daubechies Ten lectures on wavelets (SIAM, Philadelphia, 1992) 

[6] M. Farge, Annu. Rev. Fluid Mech. 24, 395 (1992) 

[7] J. Eggers and S. Grossmann, Phys. Lett. A 153, 12 (1991) 

[8] T. Vicsek and A.L. Barabasi, J. Phys. A (Math. Gen.) 24, L485 (1991) 
A.L. Barabasi, O. Szepfalusy and T. Vicsek, Physica A 178, 17 (1991) 

[9] R. Benzi, L. Biferale, A. Crisanti, G. Paladin, M. Vergassola and A. Vulpiani, Physica D 65, 
352 (1993) 

[10] A. Juneja, D.P. Lathrop, K.R. Sreenivasan and G. Stolovitzky, Phys. Rev. E 49 



[11] P. OUa, Phys. Rev. E, (in press), Los Alamos preprint chao-dyn/9707023 



[12] A. Arneodo, J.F. Muzy and S.G. Roux, J. Phys. II (France) 7, 363 (1997) 

[13] A. Arneodo, E. Bacry, S. Manneville and J.F. Muzy, Phys. Rev. Lett. 80, 708 (1998) 

[14] R. Benzi, L. Biferale and F. Toschi, Submitted to Phys. Rev. Lett.; Los Alamos preprint 
|:hao-dyn/9801007| 

[15] J.C.H Fung, J. C.R. Hunt, N.A. Malik and R.J. Perkins, J. Fluid Mech. 236, 281 (1992) 
[16] P. OUa and P. Paradisi, (1998) Preprint 



10 



